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Using Schwinger's quantum action principle, dispersion relations are obtained for neutral scalar 
mesons interacting with bi-local sources. These relations are used as the basis of a method for 
representing the effect of interactions in the Gaussian approximation to field theory, and it is argued 
that a marked inhomogeneity, in space-time dependence of the sources, forces a discrete spectrum 
on the field. The development of such a system is characterized by features commonly associated 
with chaos and self-organization (localization by domain or cell formation). The Green functions 
play the role of an iterative map in phase space. Stable systems reside at the fixed points of the 
map. The present work can be applied to self-interacting theories by choosing suitable properties 
for the sources. Rapid transport leads to a second order phase transition and anomalous dispersion. 
Finally, it is shown that there is a compact representation of the non-equilibrium dynamics in terms 
of generalized chemical potentials, or equivalently as a pseudo-gauge theory, with an imaginary 
' charge. This analogy shows, more clearly, how dissipation and entropy production are related to 

the source picture and transform a flip-flop like behaviour between two reservoirs into the Landau 
problem in a constant 'magnetic field'. A summary of conventions and formalism is provided as a 
basis for future work. 
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00 ■ I. INTRODUCTION 

o 

' An increasing number of actual problems in physics find their natural expression not in the the static (equilibrium) 
aspects of quantum systems, but in the kinematical (non-equilibrium) development of their average properties. Ex- 
i-Q ■ amplcs include studies of early universe expansion , heavy ion collisions and the postulated quark-gluon plasma 
I , lasers and other driven systems (5]-|7|] , and particle creation in changing fields of force |||| . 

The migration from equilibrium to non-equilibrium involves a shift of paradigm. In common with zero-temperature 
field theory, particle systems at equilibrium are often treated by a scattering formalism, with an initial (in) state and a 
final (out) state; this is only sensible if both are known and are at equilibrium with the same thermodynamic reservoir. 
The physics of a non-equilibrium system demands different boundary conditions. The initial and final states are (by 
definition) not characterized by the same ensemble and it is more appropriate to define the state (spectral profile or 
density matrix) of the system at some initial time U and compute the final state and its consequence at a later time 
t. This describes to an initial value problem which is deftly handled by Schwinger's closed time path (CTP) action 
principle. The new picture also implies a concern with probabilities, or expectation values rather than amplitudes. 

In equilibrium, one is used to the notion of translational invariance in space and time, implying that physical 
quantities only depend on the differences of coordinates x — x' . When the field is driven into disequilibrium, it 
acquires an additional dependence on the average position and time x = ^(x + x'). This is measured relative to an 
initial point of reference Xi. In practical applications it is usually necessary to assume that the dependence on the 
average coordinate is quasi-static or of low adiabatic order in order to make computations tractable. The dependence 
on the average coordinate has important features: the preservation of unitarity demands that the statistical state of 
the field only depend on x and not x — x' . Since the state of the field can only be altered by the intervention of 
sources or sinks (hereafter referred to collectively as sources), the sources must also develop with respect to the average 
coordinate. Since one is interested both in fluctuations and the average kinematics, it is convenient to work with 
variables and sources which are bi-local objects rather than working with the field itself. This is in contrast to the pure 
field approach used by Schwinger [fl0[ . Self-interacting theories are a special case in which the field is its own source; 
they pose mainly calculational problems — conceptually no new issues are introduced other than self-consistency. 

Since the external sources affect the average state, they can be regarded as thermodynamical reservoirs, with the 
caveat that they must suffer a 'recoil' or back-reaction as a result of their effect on the system. This is not negligeable 
off equilibrium. Many quantum systems (the laser, for example) can be treated as two-reservoir systems in which the 
'external' reservoir is of comparable magnitude to the local one. 
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The nomenclature 'open system' is used to describe a system coupled to independent sources. The name 'closed 
system' is given to a system without sources or one in which the field is its own source; in the latter case, the source 
must become effectively impotent with regard to the further development of the average state, as equilibrium is 
approached. Equilibrium is only achieved when, either the whole system approaches some driven limit cycle, or the 
contact with external sources is effectively terminated. In a closed system, the final equilibrium is a thermal state, or 
a state of maximum entropy. Some authors define equilibrium to mean a thermal state rather than merely a static 
one — and non-equilibrium to mean anything else. This is somewhat misleading since a non-thermal but static state is 
still final in the absence of new perturbations and must therefore be considered a point of equilibrium for the system. 

The purpose of the present paper is to extend Schwinger's method of analysis to treat non-equilibrium ensembles of 
bosons. This touches on and extends a number of apparently different approaches to non-equilibrium [p^Jllfl . Since 
Schwinger's original work JlO| on the initial- value problem, most authors have been seduced by the functional integral 
and have therefore missed the often subtle advantages of Schwinger's methodology. It is intended that the present 
work should convey a pedagogical flavour of the suggested approach, which overlaps with the existing literature in 
strategic places without actually following any of them. In particular, conventions and definitions (which differ from 
most other accounts) have been chosen rather carefully for practical purposes. Some well known results are rederived 
in order to make the present work as self-contained as possible. The paper begins with a summarial discussion of 
the formalism, paying special attention to the action principle and unitarity; later the most general quadratic theory 
which maintains unitarity is presented and the Green functions are calculated for prescribed sources. Particular 
attention is paid to the effects of non-locality in the sources — an issue which has been largely neglected in previous 
work, and turns out to place strong requirements on the behaviour of stable systems. Finally a brief comparison is 
made between the present work and other approaches. 



II. FORMALISM AND CONVENTIONS 



The conceptual framework for the decription of non-equilibrium processes will include operator field theory, the 
method of sources and the local momentum space Green functions. In addition it proves convenient to use Schwinger's 
quantum action principle. This is a statement about the unitary development of the field with respect to the variation 
of certain variables. Since it embodies the equations of motion and the fundamental commutation relations for the 
field, it is both compact and elegant. One begins with the action operator, which is defined to be the classical action 
with the classical field replaced by the field operator, together with a suitable ordering prescription for the fields. 
Here the ordering will be the usual time-ordering and the action that for a real scalar field without self-interactions. 
The Minkowski metric-signature is ( — !- + + •••) which allows straightforward comparison with the Euclidean theory. 

S = J dV x {±(d»<t>)(d„<t>) + imV - J4>} (1) 

where dV x is the Minkowski volume element. The operator equations of motion now follow from the quantum action 
principle jl2| 

6(<p\<f>') = i{4\6S\<t>') (2) 



giving 



(-□+m 2 )<Xx) = J(x) (3) 



Given that 4>{x) = 4>(xi) at initial time U (or, more generally, on the the space-like hypersurface Cj), the solution to 
(||) may be written 



<f>(x) = <t>{xi) + / dV x >G c (x,x')J(x') (4) 

J (7 i 

where ai and a are the initial and final hypersurfaces and G c (x,x f ) is a Green function which satisfies 

(-n+m 2 )G c (x,x') = S(x,x'). (5) 

Both the Feynman propagator and the retarded Green function have this property. 

The surface integral under the variation of the action vanishes independently implying that the generator of in- 
finitesimal unitary transformations on the field is [[L2| 
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X4> = J da^d^. (6) 

Since it is easily established that the unitary variation of any operator A is 

5A = -i[A iX A] (7) 

it follows that, on any spacelike hypersurface with orthogonal vector n M , one has 

[0,H /J ]n' 1 =^(x,x') (8) 

with = d^4>. This is the covariant statement of the canonical commutation relations for the field and its conjugate 
momentum. To avoid unnecessary notation it is convenient to write this simply as 

[0(x,i),9 t 0(x',i)] =W(x,x') (9) 

with the understanding that general covariance is easily restored by introducing a suitable time-like vector. 

From the action principle (g) it can be shown by repeated functional differentiation with respect to the source that 

' V;!i ' J " L ' / -i*{<h\T(<i>{xi)...<t>{x n ))\<h), (10) 



SJ(xi) . . . 5J(x n ) 

thus the Taylor expansion of the amplitude may be written in the shorthand form 

(tf»\<h)j = (faFeWfa) (11) 

where T denotes time-ordering (latest time to the left). This formula may be regarded as a generating functional for 
the n-point functions of the theory. The complex conjugate of this quantity is 

{<h\4i% = (4>i\<t>2)j = (^|Tt e -'^ 2 ) (12) 

where stands for anti-time-ordering (latest field to the right). This reverse-ordering is necessary to ensure the 
cancellation of intermediate fields in the identity: 

(0 2 |0i>./ x (<f>i\<h)j = (th\ r ^er iJ *Te iJ *\(h) = 1. (13) 

This is the key observation for the construction of the expectation values. Notice how the operator ordering in jl^ ) 
starts from an early time, increases to a final time (at the centre of the operator product) and then reverses back 
to the initial time. Each field, at each instant along the closed time path has a mirror counterpart required for the 
cancellation of the intermediat e op erators in (|l3|). This property can now be used to advantage to construct a 'closed 
time path' action principle Jl0|Jl3[l . 

Consider an expectation value of the form 

(t\X(t')\t)=J2(t\i)(i\X(t')\i'){i'\t) (14) 

where the sum over intermediate states i, i' is a sum over all states and (t\ is a shorthand which refers to either a 
pure state of the system, or a mixed state, specified at time t. The expectation value specifies the average value of 
the operator X at the time t' given the state of the system at time t. It involves conjugate amplitudes and hence the 
conjugate forms of the action principle: 

6(t\i)=i(t\SS ti \i) (15) 
6(i\t) = -i(i\SSj t \t) (16) 

S ah = [ dtL. (17) 



To obtain ( pT| ) from an action principle one would therefore like to introduce the operator X by functional differen- 
tiation with respect to an appropriate source (or combination of sources) between an amplitude and its conjugate. 
This is achieved in the following way. First one observes that 
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5(t\t)=dC£mx(i\t}) = 6(l) = 0, 



(18) 



so differentiation of this object is to no avail. However, if we make an artificial distinction between the amplitude and 
its conjugate by labelling all objects in the former with a + symbol and all objects in the latter with a — symbol, 

5{t\t) = lim i{t\SS+ - SSl\t) (19) 



then we can use the solution of this quantity as a generating functional for (14) since X can be expanded in terms of 
either tJ- otjj—. This breaks the symmetry of symbols in ([l8]). At the end of a variational calculation one removes 
the + and — symbols restoring the conjugate relationship between the two amplitudes, having inserted the appropriate 
operators by differentiation with respect to the source of only one of them. Note in ( |l9| ) that, for any unitary field 
theory, the action is self-adjoint, thus we may drop the dagger symbol in future. Also, in treating the + and — parts 
of the field as being artificially independent, the condition 

)=((>- (too) (20) 

is required to ensure that the limit + — > — restores the single identity of the field operators, and additionally one 
must have that all 0_ fields (at any time) must stand to the left of all 0+ fields (at any time). Since — fields are 
anti-time-ordered and + fields are time ordered, this condition arises naturally and ensures the triviality of (filf). 
The meaning of the above procedure can be illustrated by noting that the solution to ( |l9| ) may be written 

(ti\ti) J± = (t i \T^e- iJ ^-Te iJ +^+\t i ), (21) 

The expectation value of the field is found using the ordered expression 



• <5 i i \ 
oJ + [x) 



+=- 



(0(a)) (22) 



= lim (t 4 |exp ( -i j J-(j>- + i J J+4>+j </>+(ir)exp \ i J ) |t ( ). (23) 



H — 

Taking the limit H > — , one has 

(0(x)> = (t 4 |exp f-ij^ J<j\ <Xx)exp U j J A (24) 

This shows that the average value of the operator depends only on the past (retarded) history of the system beginning 
from the initial time tj. It can be shown (see appendix A) that the closed time path generating functional is closely 
related to the generator for the retarded n-point functions. The acausal (advanced) pieces cancel in the limit H > — . 

So far, the discussion has used the slightly trivial example of pure states (t\. As noted implicitly by Schwinger [[j"o|| , 
the same action principle holds when (t\ . . . \t) is replaced by (t\p(t) . . . \t) (a mixture of states) since this does not affect 
the conjugate relationship between amplitudes. The nature of the expectation value can therefore be left out of the 
discussion for the most part. Indeed, in practice, the effect of a non-trivial density matrix in the expectation value can 
be mimicked by the introduction of suitable sources [l^Jl4] | — a procedure which will be adopted in the next section. 
To present the formalism in a way conducive to generalization, the next step is to present the Green functions for the 
case of pure-state vacuum expectation values and then introduce the finite temperature (mixed state) modifications 
which will be the starting point for writing down an ansatz for non-equilibrium. 

The above use of generating functionals is closely related to the path integral approaches of Calzetta and Hu jij , and 
Lawrie S. It proves useful not to pass directly to the path integral however, but to follow Schwinger 's approach. For 
the remainder of the paper, equation ( |l9| ) will be considered the starting point for the discussion of non-equilibrium 
field theory. 

From equations ([!]) and (|l^) one obtains the operator equations of motion for the field. Taking the initial time to 
be tj, the furthermost future time to be and the final time at which expectation values are to be computed as tf, 
then using the boundary condition in equation (EOT), 



<f>+(x) = <j){x t ) + J G c (x, x')J + (x')dV x > 
<f>-{x) = Mtoo) + [ 1 G c (x,x')J-(x')dV x . 



f 



Xi )+ G c (x,x')J + {x')dV x , + G c {x,x')J-{x')dV x < (25) 
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where G c (x,x') is a retarded (causal) Green function. Notice that, as the distinction between + and — is removed, 
these equations reduce to ([|). Substituting these into the exponential solution to (|n]) and defining a vector and its 
transverse by J T — (J+ </-), one may write 



too rt»_ 

T i 



ln(0,ti\0,ti) = -i J-{x')dV x ,<f){xi)+i J T (x)G(x,x')J(x')dV x dV x , (26) 



where 



^ttM-^^^ 1 '^^ ^ (97) 

^ x > x >-{ -G c (x,x') 6(x' -x)G c (x,x') J ' V> 

and 6{x — x') is the step function which satisfies 

9(x) + 6{-x) = 1. (28) 

As a result of this property, the sum of rows and columns in ( p7| ) is zero. This is a reflection of the triviality of 
equation (|l3|). ft further implies the causality of expectation values derived from this generating functional. While 
( |27| ) has a simple physical derivation in terms of the equations of motion, a more symmetrical form can be obtained 
by attaching a variational interpretation to G c (x,x') directly. Again, following Schwinger, and varying with respect 
to the sources 

5 2 6 1 (t i \t i ) = (i) 2 (t\(6 2 S+ - 5 2 S-){6 1 S+ - fcS_)|)t > (29) 
where, according to the ordering rule, this equals 

hh{U,ti) = (i) 2 (t|<?(> + (z2)<M^i) + <P-(x 2 )4>-{x 1 ) - <Mx 2 )(Mxi) - 4>+(x 2 )<t>-(xi)\)t > . (30) 
Comparing the solution of this to 

expQ J dV x dV x <J T G(x,x')J(x')J (31) 



one has 



where 



Gi, -.,-'): ( % + _ + + % + J ). 132) 



{<j> + {x)<j> + {x')) = -iG ++ {x,x') (33) 

{<j> + {x)<j>-{x'))=iG + ^{x,x') (34) 

(cj>-(x)<f> + (x , )}=iG- + (x,x') (35) 

(<f)-{x)(f>-(x')) = -iG—(x,x'). (36) 

As the distinction between + and — is lifted, the assumed ordering implies that 

G++(x,x') = i(T<j)(x)<t)(x')) = G F (x,x') (37) 
G+_(x,x') = -i(ct>(x)(t>(x')) = -G ( -\x,x') (38) 
G- + (x,x') = -i(<f>(x')(t>(x)) = G {+) {x,x') (39) 
G—(x,x') = i{T^{x)4>{x')) = Gaf(x,x') (40) 

(41) 

where Gf is the Feynman propagator, G^ are the positive and negative frequency Wightman functions and Gaf is 
the anti-time ordered propagator. In the limit of zero source, these quantities satisfy the equations 

{-U + m 2 )G F {x,x') = 6(x,x') (42) 
(-□+m 2 )G (±) (£,£') = (43) 
{-\J + m 2 )G AF (x,x') = -S(x,x') (44) 

(45) 
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by virtue of equation (||) for the field operator. The non-zero right hand side of ( pL2f ) and (|43|) are due to the time 
ordering. From the time-ordering (|37]) it follows that 



- iG F (x, x') = i6(t - t')G {+) 0, x 1 ) - i6{t' - t)G ( - ] (x, x'). (46) 

Substituting this relation into ( |37| ) and using the commutation relations for the field (^|) proves (|37|). Similarly ( |4^ ) 
follows from the relation 

G AF (x,x') = -G F (x,x'). (47) 

Using (^8|), it is now straightforward to see that the sum of the rows and columns in ( |32| ) is vanishing, as required for 
causality. A number of additional relations between the Green functions can be proven. The retarded and advanced 
Green functions satisfy: 

G ret (x,x') = -6{t-t')[cj ) (x),<j ) {x')] 

G adv (x, x 1 ) = 6{t' - t)[<j>{x),4>{x')]. (48) 
Also, in virtue of (^) it is easy to see that 

G F = Gret + G<-> = G a dv ~ (49) 

Gap = = G ret + G<+> = G adv - G((-)). (50) 
The unequal-time commutator and anti-commutator Green-functions are defined by 

G(x, x') = [<p(x), <f>(x')] = G (+) + G ( - } 

G(x, x') = {(f>{x), <t>{x')} = G (+) - G ( -> (51) 

These will be useful later and serve to pin-point the conventions used in this work. Before considering the momentum- 
space representation of these functions it is useful to note that G{x, x') can be written entirely in terms of the formal 
quantity 

H(x,x , ) = i{t\<j>(x)(j)(x r )\t) (52) 

as 

G M = «< - ( -0$ ) * « - •) ( ) . (53) 

where 

H(x',x) =H(x,x')* (54) 

and the sum of rows and columns is manifestly zero. Since the spectrum of the operator — rj + m 2 on the complex 
wave e lkx is solved for any k satisfying a dispersion relation, the solution to ( f43| ) is the most general linear combination 
of plane waves satisfying the dispersion relation k 2 + m 2 = 0. This implies that the vacuum positive and negative 
frequency Wightman functions can be written, in n spacetime dimensions, 

/An L, 
JULjk^-xTe^Wk 2 + m 2 ) (55) 
(2lT) n 

/rl n b 
AA. e ^- X r 9{ _ ko)S{k 2 + m 2y (56) 



Defining the fourier transform of G(x,x') by 



and using the integral representation 



An i. 



f +oc duj e -<"C*-0 

9(t -t')=i — — 58 

' , 2tt u> + le 



G 



it is straightforward to show from (46) that 



G ^ = WT^Te > (59) 

which is fully consistent with (f42|). GafQs) is easily obtained from (f47j). 

Note that, had the symmetrical form of G(x, x') not been used, similar results could still have been obtained. It 
is possible, in the manner of a symmetry transformation to redefine the Wightman functions so that positive and 
negative frequencies are mixed. This simply mixes up the Feynman and anti-Feynman propagator also. For instance, 
if one defines 

G (+) (fc) = 2iri5(k 2 + m 2 )[0(k o ) + a6(k ) + f36(-k )}, (60) 
then the corresponding Feynman propagator becomes 

= * + 2 Q . - l2 J^- (61) 

where the last term is evidently a piece from Gaf- Since this only complicates the matter, such redefinitions will not 
be pursued further. 

So far, this summary of the action principle has only explicitly encompassed pure state expectation values, which 
are comparatively trivial. A statistical system with real particle densities, as well as perhaps a temperature and 
entropy, is described by a mixture of such vacuum expectation values (since the character of the actual pure state is 
not usually knowable), with the the statistical weight given by the density matrix p. The simplest example of such is 
a system in thermal equilibrium (p = exjp(—[3H)). Although a thermal system is quite extraordinary as many particle 
systems go, it serves as a useful reference point, both from the viewpoint of formalism and from a physical perspective, 
since very many physical systems can be characterized by a temperature of sorts. A statistical expectation value for 
some operator X may be written 

_ Tr(t\p(t)X(t')\t) 

(t\X{t)\t).= Tr{tlpm (62) 

and characterizes the average value of X at the time t' given the state of the system at time t. Notice that the trace is 
over probabilities of the form (t\t) rather than amplitudes it\t'). The latter would be meaningless. The structure of the 
expectation value is therefore simply that in equation (|14| ) and the closed time path action principle applies. Indeed, 
it is noteworthy that the density matrix itself is merely an operator which can effectively by introduced into the pure 
state generating functional by functional differentiation with respect to an appropriate source. There is therefore no 
loss of generality in taking the closed time path action principle at face value and making no special reference to p.. 

The cyclic property of the trace in (^2|) has noteworthy implications for the Green functions and sources in the 
CTP formalism. Consider the expectation value in (|6^). This may be rewritten as 



Tr{p(t)e iH ^'-^X(t)e- iH ^-^) 

(63) 



Tr(t\p(t)\t) 
Tfc(e- wr ( t '- t >p(t)e ar <*' -*>*(*)) 



Tv(t\p(t)\t) 

where H is the Hamiltonian of the system. Using this 'relativity' between the time-dependence of p and X, it is 
possible to place all of the dynamical development of the system in either one or the other. An example of the use of 
density matrix time-development is given in ref. [ pr| . In the CTP formalism, the distinction between forward moving 
times and backward moving times makes equation (p3|) effectively 



{X{t ))s Mt\p(t)\t) ■ (64) 

The cyclic property of the trace therefore implies that the density matrix p always sits between the + and — branches 
of the operator product and hence it must always be reflected by the off- diagonal terms in ±-space. In the special 
case of a thermal density matrix, the same observation leads to the well-known KMS condition |1J|, by identifying 
the inverse-temperature (3 with imaginary time. This is seen by considering the thermal Wightman function 
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(+) _ Tr(t\e-^0(x)0(x')\t) 

f 3 [X > X) - Tv(t\e-^ H \t) 

_ Tr(^|e-^ g 0(x)e^e-^^(a: / ) |t) 

~ Tr{t\e-P H \t) 

_ Ti(t\e-^ H (^(x')^(x,t + i0)\t) 

~ Tr(t\e-^ H \t) 

= -G { -\n,t + il3,x') (65) 

where the cyclic property of the trace has been used. The left and right hand sides are precisely the elements of the 

off-diagonal G-\ and G |_. This property is, in fact, sufficient to determine the thermal Green functions. 

To determine these, in a form which manifestly reduced to the vacuum case, one writes 

G (+) (fc) = -2m[9{k Q ) +X}S{k 2 + m 2 ) 

G (_) (A;) = 2m[6(-k ) + Y]5{k 2 + m 2 ) (66) 

with X and Y to be determined. Since the commutator function G(x, x') must be independent of the state of the 
system (in order to preserve the canonical commutation relations), it follows immediately that X = Y. If one then 
employs the KMS condition which, in momentum space, becomes the definite relation 

G {+ \k) = -e^G^ik) (67) 

it follows that 

9(ko) + X = e^[9(-k ) + X] (68) 

and hence 

X = 9(k )f(\k Q \) = f >a (k ) (69) 

where 

By considering (amongst other things) G(x,x'), it follows that /(|feo|) is an even function of fco thus f(ko)9(ko) = 
f(—ko)9(—k()), whereupon it is trivial to show the unitarity relation 

G {+ \x,x') = G ( ^(x,x')*. (71) 

Note that the fact that G^ consists only of positive frequencies fco is pivotal in this derivation. The Feynman 
propagator can now be obtained from equation (|4^) by using the integral representation of the step- function (|5§|). 
The thermal Green functions are therefore summarized by 

G (+) (fc) = -2m6{k ){l + f{\k \)]6(k 2 + m 2 ) (72) 
°F( k ) = ,2 i 1 2 + 2mf(\k \)S(k 2 + m 2 )9(k ). (73) 

Another important form of G^ is obtained by performing the integral over fco, thereby enforcing the role of the 
delta-function in ( |66j ) . This gives a result which will prove more useful for calculations later and is more closely related 
to the ansatz used by Lawrie in ref. ^ : 

G^(x,x') = -2ni ( /^ e *(Mx-x')-^-0) + 2 />o(M)) 
J (27T)™- 1 2\w\ 

--27TJ [ dn lk r '(k-(x-x')-^(t-t')) ( 1 + /(l a; l)) f741 

2m J (27T)-! 6 2| W | 

where />o(fco) is the function composed of only positive frequencies. It is now straightforward to verify that the 
canonical commutation relations are satisfied, by differentiating G(x, x') with respect to t' (see equation (|5l])). 

This completes the presentation of conventions to be used in the remainder of the paper. It is convenient to add 
here that a bar (e.g. a) represents an object which is even, while an object with a tilde (e.g. a) represents one which 
is odd with respect to its arguments. 
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III. INTERACTION WITH SOURCES 



The formalism demonstrated so far has been for free fields. Free fields are always in a state of equilibrium and 
therefore the discussion needs to be widened to incorporate collisions or interactions. The present work will deal 
with interactions which can be mediated by sources of the type (j>(x)A(x,x')<j)(x'). This includes a variety of self 
interactions, contact with external forces and noise or impurity scattering, depending on the nature of A(x,x'). The 
self-energy of an interacting field theory has this form, for instance, thus sources of the quadratic type can also be a 
representation of the lowest order, self-consistent 'particle dressings'. In Feynman diagram language, these represent 
the re-summed one-loop approximation, Hartrce approximation and so on. Lawrie || uses the notion of such sources 
to effect a renormalization (resummation) of perturbation theory in a real scalar field theory. The same idea is 
expressed in a different language in the work of Calzetta and Hu |l[ . Since it is not the aim of this paper to discuss 
specific models, the specific nature of the source terms will not be specified here. Rather the discussion will centre 
around what general properties such a system might have and a discussion of possible applications will follow. 

In an interacting theory, one normally perturbs about the free field theory. Unfortunately, the dispersion relation 
(or 'mass shell' constraint) for free particles is no longer appropriate, since it reflects none of the interactions which 
'dress' the particles. A more satisfactory starting point would be a 'quasi-particle' mass shell, including some of the 
interaction effects as the basis for a perturbation theory. This is the essence of a renormalization and can be effected 
by the use of sources jl6| . 

The starting point for the investigation of non-equilibrium fields will therefore be the closed time path (CTP) action 
principle, taking the action S for a neutral scalar boson and supplementing it by quadratic sources. Observing the 
CTP operator ordering, one has 

S+ - S- - S+ - S- + I [ dV x dV xl [T{4>+{x)A ++ {x,x')4> + {x')) 



2 t 

+ <j>-(x)A_ + (x,x')(f> + (x') 
+ <t)-{x')A + ^(x 1 x')4) + {x) 

+ T\d>.(x)A^(x,x')^(x'))} (75) 

It should be clear that no fundamental field theory may contain off-diagonal terms in ±-space. The CTP action 
principle is, by construction, diagonal, being the difference between S+ and 5_ (see equation (|l9|) ) . However, it 
was remarked earlier that the effect of a density matrix must be reflected in off-diagonal terms so, while such terms 
are certainly not fundamental, they can exist as off-diagonal self-energies representing the dynamics of a density 
matrix. Moreover, since off-diagonal terms represent a point of contact between fields moving forward in time and 
fields moving backward in time, one might anticipate that off-diagonal sources would be at least partly responsible 
for choosing an arrow of time (the generation of entropy) . The explicit coupling will therefore play an important role 
in both non-equilibrium kinematics and dynamics. 



The essential unitarity of the CTP formalism is seen, from equations (15) and (jlq), to be summarized by the 
following property of the transformation function: 

(ti\ti)± = (tilhh ( 76 ) 

namely that complex conjugation merely exchanges + labels with — labels and vice versa. If one defines indices 
a, b = +, — , then the operator defined by the second variation of ([75|) with respect to the field <p a , S a t = S a 5b(S+ —S-), 
with S a b = S++, S-\ . . ., satisfies the relations 

Si_ + (x, x') = —S--(x,x') 

S^_(x,x') = -S-+(x,x'). (77) 
This, in turn, implies that S ao may be written in terms of real constants A, B, C and 7^. 

s , x ,s = 1 (-□ + m 2 )5{x, x 1 ) + A(x, x 1 ) + iC(x, x>) B(x, x') + ^(x, x 1 ) ffr -iC(x, x') , ( ^ , 

-B(x,x') -Y{x,x') -iC{x,x') qz\-m 2 )8{x.x') -A(x,x') +iC(x,x') 

where a new derivative has been defined to commute with the function 7^(x, x'): 

Dl = dn +^lv 1 (x,x') 9 M J v (x,x'). (79) 
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Note, first of all, that the sum of rows and columns in this operator is zero, as required for unitarity and subsequent 
causality. Derivatives higher than first order in the sources could be rewritten using the field equations (still to be 
found) and absorbed into other terms, thus such terms are redundant. There can be no non- vanishing terms of 
the form + i90 + without violating time reversal invariance or merely adding total derivatives to the action. Finally 
C =/= is clearly disallowed in a fundamental theory on the grounds of unitarity. It turns out, by considering the 
field equations, that the only fully consistent choice is C = 0, even though such a term does not violate equation 
(f76|). Equation ( |78| ) agrees with the form given by Lawrie, up to differences in conventions and the inclusion here of 
B(x,x'). 

The significance of the off-diagonal terms involving 7 M can be seen by writing out the coupling fully: 

7 M (*,Z')- (01^2-02^01). (80) 

The term in parentheses has the form of a current between components 0i (the the forward moving field) and 02 
(the backward moving field). When these two are in equilibrium there will be no dissipation to the external reservoir 
and these off-diagonal terms will vanish. This indicates that these off-diagonal components (which are related to off- 
diagonal density matrix elements, as noted earlier) can be understood as the mediators of a detailed balance condition 
for the field. A similar conclusion was reached in reference by a different argument for the quantity referred to 
here as B. When the term is non-vanishing, it represents a current flowing in one particular direction, pointing out 
the arrow of time for either positive or negative frequencies. The current is a 'canonical current' and is clearly related 
to the fundamental commutator for the scalar field in the limit H > — . 



Using equation (78) it is now possible to express ( |75| ) in the form 

S C tp = j dV x dV x ,^ a S ab dp b (81) 
and thus the closed time path field equations may be found by varying this action with respect to the + and — fields: 
^gg = (_ D + m 2 )0+(a . )+ l J dK ,|(3 + i c) 0+ ( x ' ) + (B_,c)0_ (x ' )+ ^(a^_( s O)+^0-(x')| (82) 

SSctp 
5(j>-(x) 



= O-m 2 )0_(a;)-- J dV x , I (A - iC)<t> + {x') + (B + iC)<f>. (x') + Y{d^+{x')) + c^ 7/i + (x')j (83) 

(84) 

and setting the right hand side to zero, which introduces the notation 

A(x, x 1 ) = - (A(x, x 1 ) + A(x', x)) 

A(x, x') = - (A(x, x') - A(x', x)) 

1 

2 



d^(x,x') = - ( 9" Mx,x')+ d» 7fl (x',x) ) . (85) 



The canonical commutation relations which derive from the action principle (see equation (^)) are unchanged by 
these modifications, since they depend only on x and therefore cancel in the commutator. This is key feature in any 
consistent description of non-equilibrium phenomena. 

To solve this system of non-local equations, the best strategy will be to look for the Green functions, or the inverse 
of the operator S a b- This is the method adopted by Lawrie ||. Although a common strategy will be used here, the 
method will be somewhat different in spirit. The variational approach used in rcf. jl| will not be used here. Owing 
to the non-locality, it is clear that the inverse operator cannot be a translationally invariant function in the general 
case. It must be formally dependent on both Cartesian coordinate differences and the average coordinate: 

x = ^(x-x') (86) 

x=l(x + x'). (87) 

Moreover, since the operator contains off-diagonal terms, which typically signify a non-trivial density matrix, it is 
natural to look for a solution based on the form of equation (f72h , generalized to include a dependence on x. Although 
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this sounds like an innocent step, it is far from a trivial undertaking since it introduces non-linearities in the spectrum 
of excitations which must be handled in a self-consistent way. It is useful to work with the quantity H{x,x'), from 
which all the Green functions can be obtained (actually the Wightman function in disguise). Using either the field 
equations or the matrix equation 



S ab G bc = 5 a c S(x,x') 

one obtains equations of motion for the quantity H(x,x') (see equation (p3|)). Not all of these equations are inde- 
pendent, owing to the symmetry in equation (uw. In particular, their consistency requires that C = which is now 
chosen explicitly. It is sufficient to consider 

{-U + m 2 )H{x,x') + j ' dV x n h{x,x") - B(x,x') -^(x,x") % » -d^(x,x")j H(x" ,x) = (89) 

on the understanding that H(x,x') — H(x',x)*. This relation is to be supplemented by the canonical commutation 
relations for the field, which appear in equation (p8|) in the form 



d t > (H(x,x')-H(x',x)) 



iS(x,x') (90) 



and complete the consistency of the equations of motion. 

The next step in obtaining an intuitive formalism is to introduce a (local) momentum space technique by Fourier 
transforming x and retaining a dependence on the average coordinate x: 

-±±e ik ^-^H(k,x). (91) 

A suitable ansatz for this function, which generalizes the dispersion relation and the one-particle distribution function 
f(ko) is 

H{k,x) = 2n0(ko)[l + f(k Q ,x)}6(-k 2 + oj 2 (k,x)). (92) 

The spacetime dependent function f(k,x) is often referred to as the Wigner function and signifies the inhomogeneity 
in particle occupation numbers. The generalized dispersion relation takes generic form — k 2 . + lu 2 = 0. In the free 
particle limit uj 2 = k 2 + m 2 . It is the determination of this dispersion relation which is of specific importance, since 
this determines the spectrum of excitations for particles in the plasma-field, and forms the basis of all perturbation 
theory when the sources represent self-interactions. 

It can be verified that, since H{k,x) depends only on the average coordinate, the commutation relations are 
preserved (see equation (|o|)). A more useful form of (^) is obtained on performing the integration over fc . This 
eliminates the dubious derivative of the delta-function from subsequent relations and leads to a number of helpful 
insights. 

H {x, x') = 2n[ e tt *<— + (93) 

v ' J (2n) n - 1 2\u\ y ' 

where it is understood that ko = \u>\. Finally, it is useful to define the derivative with respect to the average coordinate 
d = | (d x + d x > ) and the quantities 



F„ =YTf = l F ^=l^M + II (94, 



d uj 1 1 

n„ = = = In (95) 

LU Z L 



IV. DISPERSION RELATIONS 



To determine the dispersion relations for given sources it is useful to distinguish three cases which will be referred 
to as the local, translationally invariant and inhomogeneous cases respectively. In the local case, the sources are 
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proportional to a delta function. In the translationally invariant case A(x,x') = A(x — x') and in the inhomogeneous 
case A(x,x') = A(x,x). 

There are two ways in which one can proceed with the determination of the dispersion relations. One is to separate 
real and imaginary parts and the other is to used complexified momenta. The latter has several advantages and 
makes straightforward contact with the classical theory of normal modes. It will be used exclusively for determining 
the spectral relations. Separating real and imaginary parts on the other hand allows one to identify imaginary 
contributions as a Boltzmann/Vlasov equation, illustrating nicely the intimate relationship between transport and 
dissipation jL4|. 

In order to extract information from the equations it is necessary to undertake an approximation scheme in which 
only low order derivatives are kept in x. This is equivalent to an adiabatic (or quasi-static) scheme in which the 
development of the system if slow in comparison to fluctuations, so that fast and slow moving variables separate 
in an assumed way. In fact this is already built into the assumed form of the solution for the Green function, since 
without such an assumption, there are no grounds for assuming that x and x would separate in the prescribed manner. 
For most purposes this approximation is quite sound. For the present, there seems to be no way of eliminating the 
approximation. 



A. Local sources 



In the local case, the equation satisfied by H(x,x') is 



[-□ + m z + A(x) - - 7^] = 0. (96) 

Note that, since B is an odd quantity, it does not appear in the local limit. Since one is interested in the variables 
x — x' and x, it is convenient to Taylor expand x around x. Under the Fourier transform this takes the appearance 

A(x)H(x,x') - [A(x) + kd^A)^- + . . .]H(k,x). (97) 



It is useful to define a new quantity by 



2 V " ; dk„ 



df 

T> = ^ (98) 

(the steepness of the spectral envelope for the Wigner function) so that 

1 c) M 

where v g = ^ is the group velocity of the dispersing wave-packets. In terms of the quantities (Q) and (jo^) the action 
of the spacetime derivative operator on H(x,x') gives 

d, H(x,x') = 2irj JL-JLjK*-* + F„ - (100) 

and subsequent derivatives are obtained in a straightforward way. 

Substituting H(x,x') into the equation of motion ( |96"| ) now leads to a second order differential equation for the 
frequency lj{x): 

u? - 2{ik» + F» + l-T)^ + n 2 =k 2 + m 2 +A(x) + J(S M 3)(T" - </w) 
2 2 M 



- c>^ -F- IiWF^ - d^F^ - zA:% - (101) 

Clearly this equation presents an insurmountable problem for the purposes of analytic calculation, thus an approx- 
imation must be made, based on the adiabatic evolution of the average properties of the system. The lowest order 
case (which will be sufficient to reveal the features of interest in this paper) is when f2 M and derivatives of F^ may be 
neglected. This corresponds to a near classical transport of particles, with relatively few of the fluctuations added by 
the quantum nature of the field. With this approximation the dispersion relation may be written: 
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k 2 + m 2 + A(x) + -(0 M A)(T" - v^/uj) - — F 2 — 2ik^ - i^k^ - = 0- (102) 

This can be separated into a more appealing form as the dispersion relation for a damped oscillator array 

- oj 2 - iTuj + uj 2 = R (103) 

where one identifies the natural frequency, 

uj 2 = k 2 + m 2 - F 2 (104) 



the decay constant, 



and force term 



r = -^j(fl M A)(T" - + h»(F, i 7M ) ( 105 ) 



R = ^F^ + d^ tl - A(x). (106) 

One notices how the effective mass of the theory is reduced by the gradient of the Wigner function F^, indicating 
that rapid transport could lead to a second order phase transition. This might also lead to anomalous disp ersion. 



In a true linear oscillator array R, T and ujq would all be independent of the frequency uj. In equation (105) only the 
zeroth component of the last term is independent of uj. This indicates that the decay/amplification of certain modes 
in time is oscillator-like, but that the spatial modes are multiplied by a factor of k/w, the inverse phase velocity, 
which has a critical value when m/k is a maximum. This signifies the effect which a gap in the frequency spectrum 
can have in leading to anomalous dispersion in the 'plasma'. At high frequencies V — > k°(F — ^7o) and the system is 
oscillator-like. At low frequencies, damping is dominated by the external potential A and by transport as one might 
expect. 

B. Translationally-invariant sources 

In the translationally invariant case, all variables are a function of x — x' . One may therefore fully Fourier transform 
the sources: 

/rl n h 
f rl n h 

^(*-*') = J -^^ {x - xr T{k) 

Jt^_ e ik^x-x TiB(k). (107) 

(27T)™ 

Note that, since B is an antisymmetric function a factor of i is introduced to keep B(k) real. The equation of motion 
for H{x,x') is now 



[-n+m ]H(x,x') + J dV x »{A(x-x")-B(x-x')-Y{x-x') d „ -(d^)]H(x", x') = 0. (108) 

The translational invariance enables the latter spacetime integral to be performed immediately, yielding the dispersion 
relation 

k 2 + m 2 + A{k) - 2ik^^ - iB(k) = 0. (109) 

An apparent consequence of the translational invariance is that F^ — owing to the steady state nature of the system. 
Comparing the dispersion relation to equation (|103|), one identifies 



r = -k^ + B(k) 

UJ 



2 i,2 , 2 

uj n = k + m 





R = -A(k) (110) 
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Although the translationally invariant theory describes only steady state disequilibria, it is nevertheless seen that the 
field oscillations are concentrated around the usual mass-shell u>q with an amplitude driven by the external force 

m (HI) 



[(W 2 - W 2 ) 2 + ( rw )2]* 

and a quality factor Q = ujq/T. Such a steady-state description would be appropriate for an 'infinite laser' i.e. a 
device which is not affected by any finite size considerations. 

C. Inhomogeneous sources 

The main case of interest is when the sources and Green functions have a residual dependence on the average 
position and time. This includes the local limit as a special case: 

A(x, x') = a(x — x')(3(x + x') 
a — ► 8(x — x') 

13 ->A(x) =A(x'). (112) 

As usual, one is looking for the eigenspectrum of the quadratic operator acting on H(x,x'). The equation satisfied 
by H (x, x 1 ) is now: 



[-\3+m 2 }H(x,x')+ J dV x , 



A(x, x") - B(x, x") - Y(x, x") -8^ 



H(x",x')=0 (113) 



In the inhomogeneous case there is no dispersion relation consisting of continuous frequencies in general so the 
dispersion relation will only exist for a discrete set. It is convenient to divide the discussion into two parts: the 
determination of the dispersion relation and the nature of the restricted set of values which satisfy the dispersion 
relation. 

The problem to be addressed is contained in following form in momentum space: 

/d n h d n n 
dVx " ^e^'^+^'-^Sip, x" + x')H(k, x + x') = XH(x, x'). (114) 

The integral over x" is no longer a known quantity in general, but it is possible to extract an overall Fourier transform 
by shifting the momentum p — > p + k and defining the average variable of interest x = ^(x + x'): 

(k 2 + iWd^ + m 2 )H(k, x) + 

/d n n 
dV x n—4-S(k, x + x")H(k + p, x" + x 1 ) = XH(k, x). (115) 
(2n) n 

In order to find eigenvalues, it is necessary to extract the factor of H(k,x) from this expression. This is not possible 
for arbitrary values of k. It is possible, however, if the momenta are restricted to a denumerable set expressed by the 
property 

H(k+p,x" + x') =H(k,x" + x') (116) 

which implies that H(k,x) is a periodic function of the momenta. The absence of eigenvalues or the failure of this 
property leads to the consideration of an infinite iterative mapping of states, which — in the absence of a stable limit — 
suggests chaotic excitations of the field. This can also be argued geometrically (see the final section). Given this 
mitigating condition, one has 



/ 



d n r> 

dV x „—^-S(k, x + x")H(k + p, x" + x') = S(k, x)H(k, x). (117) 

(ZTT)" 



The dispersion relation is now obtained in a straightforward fashion, adopting the same adiabatic approximation as 
before, and is given by the implicit relation 
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k 2 + m 2 + A{k,x) + i(0 M A)(T" - v<i/uj) - iB(k) 
2 



- d^(k, x) — (F — N) - 2ik»(F - A% - 2*fk» - ^{F - A% - 0. (118) 

where it is noted that {fc} is now discontinuous. Note that the antisymmetry of B(k) makes it independent of x. 
Comparing to the oscillator equation, one has 

T = -^(5,A)(T» - v*/w) + hfiFn - + 7m ) + B(k) 
uol = k 2 + m 2 - (F - N) 2 

R = TF^+W^(k,x)-A{k,x) (119) 

where will be defined presently. 

We now turn to the consequences of the condition in equation ( 116 ). There are various precedents for such a relation: 
one is Green functions defined on a torus (finite temperature, Matsubara formalism, electron band structure); another 
is the case of Landau levels on a torus. The periodicity is clearly the important factor here. In most of these cases 
the periodicity is one in real space and the result is a discrete spectrum of eigenvalues. Here the periodicity lies in the 
momentum itself. In fact the two notions are closely related and a periodic system in real space has Green functions 
which are periodic in momentum space owing to an infinite summat ion over discrete frequencies (which is therefore 



invariant under shifts by a whole number of periods). The relation ( 116 ) must be satisfied for all legal values of the 
momentum, thus the implication is that the system is degenerate — i.e. there exist bands of energy which leave the 
Green function invariant under certain shifts. These need not all refer to the same band. It is therefore possible to 
write 

H(k) = H (^-^j (120) 

where P M is the momentum periodicity length (which has dimensions of inverse space-length). This finite length must 
diverge to infinity when the inhomogeneities vanish. There is only one natural momentum/length scale which has 
these properties, namely 

L fl = P' 1 = d fl H(k,x). (121) 



In deriving (117) we have used the fact that 

d n k 



exp(ikx) = 6{x) (122) 



Since k is now restricted to a discrete set, the correctness of this relation could now be an issue. It can easily be 
verified using the formulae 



sin(fca;) = sin 



k=l 



(n + 1) . (nx) x 

sin cosec— 

2 2 2 



E. (n + 1) (nx) x 
cos(fcx) = sin — cos - cosec— (123) 

k=Q 

that an extra finite imaginary contribution can arise from the discrete nature of the spectrum, which vanishes in the 
continuous limit. It will be assumed that such a contribution can be absorbed by a redefinition of the sources. 

Although one is looking at periodic functions, the solution for H (x, x') need not be sinusoidal. In the case of 
Landau levels on the torus |T^|l8| periodicity is only achieved at the expense of a flux-quantization condition which, 
again, involves a degeneracy of solutions. There is, in fact, an analogy to this situation here. The extraordinary 
properties of Landau levels on a torus can be attributed to the non-translational invariance of the electro-magnetic 
vector potential. The similarity here is the non-translational invariance of the many-particle state as expressed by 
the dependence on x. This point will be discussed at greater length in the final section, to avoid its meaning being 
lost in the present analysis. 

The extra terms containing can now be explained. They arise from the x dependence of the momentum space 
measure: 
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giving a contribution 

N, = y°--- £ "-i> (125) 

which compounds the non-linearity. The above restrictions have no special consequences for the Feynman propagator, 
since the nature of the momentum is not used to obtain it. This is gratifying since the Feynman propagator must 
always be the literal inverse of the quadratic part of the + time-ordered action. Only the nature of the singularity is 
altered in accord with the modified dispersion: 

G ^( fc ) = u2 , 1 2 + 2m f( k > z) 5 (- fc o + w 2 )0(fc o ). (126) 

The appearance of a natural length scale, connected to the inhomogeneities of a non-equilibrium system, is an 
important feature for two reasons. Firstly, the spontaneous generation of a length scale implies the possibility of 
domain formation, or a cellular localization in the field. Secondly, the dependence of the Green functions on themselves 
implies that the stable solutions of the system can be regarded as fixed points of an iterative map. Such maps have 
been studied in connection with classical chaotic systems p9| . 

In the present case, the function H{k,x) depends not merely on itself but on its derivative. For exponential-like 
solutions one could expect that this would amount to the same thing, up to a constant multiplier. The situation would 
then be something akin to H — H(XH), for some constant A. This bears a noteworthy similarity to Feigenbaum's 
functional equation which can be written 



g(x) = ag(g(x/a)) (127) 



subject to a boundary value, or rewriting: 



g(g(Xx)) = Xg(x). (128) 

This equation has an analytic solution as a power series 

g (z) = l + ciz 2 + c 2 z 4 + ... (129) 

where a limiting value is approached through a geometric progression with Feigenbaum ration T = 4.66 and universal 
scaling factor a — —2.5. Solutions to this equation which fall outside the fixed point behaviour can be expected to 
lead to chaotic behaviour. This strongly suggests that the non-equilibrium Green functions must exhibit universal 
behaviour or chaos in their approach to stable behaviour. In other words, the approach to equilibrium need not be of 
the simple damped or over-damped form of a linear oscillator array — it could easily entail a chaotic attractor. 

V. ENTROPY, TEMPERATURE AND THE KMS CONDITION 

For systems close enough to a thermal state, it is possible to define an approximate temperature and entropy. The 
entropy of the system may be defined in various ways, often based only on combinatorial considerations of the micro- 
canonical picture. Here it is convenient to define an 'oscillator effective entropy' which is easily related to quantities 
which arise in the analysis. Suppose the Wigner function is given by the approximate equilibrium form 

f(ko, x) = (exp(/?(aSMx)) - 1) 1 (130) 

then one has 



F, = -f 



(131) 



and, classically, the statistical entropy S is 

S = k(lnZ + (3(cj)). (132) 
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For a harmonic oscillator, one has (see for example ref. p0|]) 



lnZ=-Tr(l- e -^)-i/?<c>, (133) 
thus the oscillator entropy may be defined as 

S= ~/3<u;>+Trm(l + /). (134) 
This motivates the definition of a simple measure of entropy for the oscillator array, given by 

/d n k 
—e(k )\n(l + f)S(-k 2 +u; 2 ). (135) 

The rate of change of this entropy is then 

duS E = d, 



d n ~ x k ln(l + /) 



M J (27T)"- 1 2\u 

d n ^k (1 + /) 
(27T)"- 1 2\w\ 



k>0 

~F n - H M ] (136) 



This quantity can be compared to (100). It shows that the entropy gradient can be thought of as a 'connection' for 
the field modes. The generation of entropy is therefore fundamentally connected with the flow of particle occupation 
numbers and the 'downgrading' of the frequency spectrum — i.e. the rate at which energy becomes unavailable to do 
work. 

As mentioned earlier, the effect of a non-trivial density matrix, cither at the initial time or later, is reflected in 
the off-diagonal sources and Green functions. If one imagines that the sources A±± arise from a coupling to another 
oscillator system |T(| or that they represent the self- interaction of the field to order 4 , then A±± is essentially the 
Green function for the field concerned and one would therefore expect the KMS condition to hold for the sources at 
equilibrium — now in the form 

9(\u}\)A+-(lj) = e f, "0(-\u)\)A-+(u). (137) 

This condition does not hold in general, but for an isoentropic process, in terms of the defined quantities at 7 M = 0, 
one therefore has 

B(u) = -e 0u B(v) (138) 



It is verified that 



9{lu)A+^(uj) _ gM 
9(-u))A-+(u)) 



(139) 



giving A_| = sirih(~/3|a;|)a(a;) for some a(ui) or 



= 1+ - p/?M = tanh(-/?M) (140) 

which agrees with Schwinger's result jnl. Note that the initial state f(x~i) and its subsequent development enter only 
as boundary conditions to the Green functions and the Wigner function. The changing form of f(x) is determined 
solely by the sources A±±. Thus, if the sources do not evolve, neither does fix) and neither does the implicit density 
matrix. In the perturbation around free field theory Q, f(x) always represents the state of the system at the initial 
time. 

In the approach to equilibrium one normally expects that dependence on the average coordinate x to disappear. 
This is an expression of what is often called 'loss of memory' of the initial state, since x is measured relative to the 
initial time. An equilibrium state (thermal or otherwise) is, by its nature, either static or periodic, thus the resulting 
Wigner function f(k ,x) must either be independent of x or a periodic function of this variable. One of the advantages 
of the present formulation is that one sees how the sources are responsible for this loss of memory. Since the sources 
drive the system, f(ko,x) can never become x independent as long as the sources are x dependent. Thus equilibrium 
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will only be secured by accounting for the back-reaction of the sources to the behaviour of the field. Explicit equations 
of motion for the sources have not been considered here. 

An example of a periodic 'equilibrium' is the case of Rabi oscillations in the laser (see ref. [pi] for a review), in 
which the source and the field enter into a pendulum-like flip-flop behaviour. An example of this will be given in the 
final section. 

The decay of field modes is exponential, per mode and is mediated by the source 7**(ar, x') and the gradient of the 
potential A(x, x'). This does not preclude other behaviour for the Wigner function. For example, in the simplest case 
close to equilibrium in which the system is quasi-static and A = B = 0, with almost no external force (see equation 
( |119| )), one has 7 M ~ and thus d^F^ + F 2 ~ giving F M ~ x~ x — a 'long tail' power law decay which parallels the 
decay of harmonic waves in curved spacetime [ p2[ . 



VI. CALCULATION OF EXPECTATION VALUES 



The closed time path formalism codifies the causal relationship between source and response, for the computation 
of expectation values in a general mixed state. Since it is redundant except as a calculational aid, it's introduction 
should be justified by an example. The causality of the method is not affected by the introduction of the sources 
A±±, but the dissipative dynamics are. Normally a fundamental Gaussian theory can never show dissipation, but in 
the present situation one has sources which can drive the field modes and redistribute energy. 

There are two cases of interest. In a self-interacting theory one might identify A±± with the correlation function 
for the field itself XG±±, giving rise to dispersion relation of the approximate form 

k 2 + m 2 + \Tr(k 2 + m 2 )- 1 = 0. (141) 

This is like the variational method used in ref. fl . Lawrie H takes the view that the sources can effect a renormalization 
of a self-interacting theory by choosing them in such a way as to 'minimize' the effect of higher order perturbative 
contributions. In either case, the effective 'resummation' induced by the sources makes it possible to see damping of 
field modes at the one-loop (Gaussian) level. 

Consider the response of the field to the source J(x), in the presence of A±±. One is interested in the causal 
expectation value of the field at time t, given the state of the system at the initial time. The time dependence, in 
the present formalism is now contained entirely within the sources — or equivalently the dispersion relation. That the 
CTP generator leads to a causal result is easily verified by realizing that the expectation value of the field is always 
coupled to the sources by the retarded n point functions. For an arbitrary action 5 [</>], 



(0(x)> = -i 



SJ+(x) 



(o o> = 



= i/" dV x ,[2G++{x,x') +G+-[x,x') + G. + (x,x')]J(x') + . . . 
= i J dV x > [2G++{x, x') + G ( +) - G^]J + ... 
= [ dV x ,[G ++ (x,x')-G^}J{x') + ... 

dV x >G ret (x,x')J(x') + ... (142) 



thus the expectation value depends only on retarded times. Furthermore, the result is real (being a probability) since 
the retarded Green function is explicitly the real part of the Wightman functions, restricted to retarded times by a 
step function: 

G ret (x,x') = -6(t-t')[G {+) (x,x r ) + G {+) *(x,x')}. (143) 
Making use of the integral representation (|58|), one has 

/ ^^ ° p( '" ( ',;?:; k(x ' x,)) ( G ' ti ' t » +G '-'w)^ ^ 

Relabelling and inserting the momentum-space forms for the Wightman functions from ([72]), one has 
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G ret (k) = - [ dw- { — S(u!-U)+) —5(w + W-) ) (145) 

J k - cj + it \uj + uj- J 

where lj + and — cj_ are the positive and negative frequency solutions to the appropriate dispersion relation. These are 
complex numbers in general, owing to the non-vanishing imaginary part labelled as Y, Now, since unitarity demands 
that (x, x') be the complex conjugate of G^~\x,x'), it is clear that 

LU* + =U!-. (146) 

It is assumed here that the dispersion relation has two complex roots. The quantity appearing in the delta function in 
equation (|9^) is then —k 2 +07+tJ_ which may also be written — k$ +u>*uj. To avoid confusion with previous notation 
for the absolute value, the complex modulus will not be denoted |cj|. This indicates that, in spite of the complex 
momenta in the dispersion relation, whose role it is to capture dissipation and transport/kinetic effects, the 'mass 
shell' constraint is real. The simplest expression for the retarded Green function is therefore 

\ZLU + {Kq — UJ + + It) 2cJ_(fco — C<J_ + It) J 

This expression is not manifestly real, since it is a momentum space result. However, if one defines 2iu) = u>- — oj + 
and 2uJ = lj + + w_ , where Cj and To are real, then it is possible to write 

1 / (iCjko — uj*ijJ + 2ZU 2 )(— fcp + u)*lq — Aikt^Cj) 
{-k% + u*w) 2 + Wk%w' 



This may be compared to equation (111) and reduces to 

1 



(149) 



when lo* = lu and e — > 0. Since the imaginary part of ( |148[ ) is odd with respect to the momentum variable k, 
the Fourier transform back to configuration space is real, as expected. The desired expectation value is therefore 
manifestly real and causal, and the time dependence since the initial time is contained entirely in the x dependence 
of the frequency to. 

VII. REFORMULATION 

In the preceding sections, it has been shown how dissipation and amplification of spectral modes can be incorporated 
into the dispersion of a quadratic theory, for suitably adiabatic processes. It is now practical interest to show 
that the same results can be presented in another significant form by introducing a 'covariant derivative' £)„ which 
commutes with the average development of the field state. This description parallels the structure of a gauge theory 
(in momentum space) with an imaginary charge. Alternatively one may speak of a generalized chemical potential for 
the 'gauge' field. 

Consider the derivative 

D„ =d li -a ll (150) 

and its square 

D 2 =□- <9^a M . - 2a^<9 M +a f "a fi . (151) 
Without any approximation, it is straightforward to show that, in the general inhomogeneous case, 

{ - {ik^ + f^-q^- n^) 2 - d^(ik^ + F fl -n tl - i\g }. (152) 

It is then natural to make the identification 

dp = - tip - + 7 M 

= d^S E (k)-N^+j^ (153) 
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where the meaning of this notation is such that the expression only defined when all objects are under the momentum 
integration — this is to be understood in all future expressions. The field a M is clearly related to the rate of increase 
of the entropy Se, the damping factor 7^ and the rarefaction of the localized cells _ZV M . One now has: 

f rl n ~ l k n -I- f\ 

(-D 2 + m 2 )H(x, x') = 2n / ,„ " L + f2 {-(ifc p - 7 J 2 ) - - 7 J} 



= 2tt 



(27T)"- 1 2|w| 1 y 
(2tt)"-i 2M ^ + 2ifc ^ ~ ^ ~ + (1M) 



Adding the appropriate source combinations for the inhomogeneous case one has, without approximation, the differ- 
ential equation satisfied by H(x,x'): 



-D 2 + m 2 + 7 2 (fc,x) + A(k,x) - B{k) + -{8^A){T^ - v%/u) 



H(x,x')=Q (155) 

k 



where the appearence of the subscript k to the bracket serves to remind that the equation only exists under the 
momentum integral. The local limit is simply 

[-D 2 + m 2 - 7 2 (a;) + A(x)] k H(x, x') = 0. (156) 

The 'gauge' field a M couples via an imaginary unit-charge plays the role of a generalized chemical potential on the 
manifold of positive energy solutions for the real scalar field (the chemical potential has no meaning for the full field, 
since particle numbers are not conserved) . Suppose now that one defines the analogue of the field strength tensor 

f^ = d^ l a v -d v a ll . (157) 

In many cases one will have / M „ = 0, thus one can 'gauge transform' the field, which maps 

4>{k) -» cj}(k)J ' a » dx " 

= ct>(k)e-f^- n »- N »-^ )dx " 

= cj){k)e~ SE -I {N - + ~^ dx \ (158) 

This shows the explicit decay (amplification) of the fc-th field mode. The latter relation shows that this process 
involves an increase in the effective oscillator entropy of the system. 

In terms of the above formulation, the spectral content of the bosonic theory reduces to the problem of finding the 
eigenvalues of the operator D 2 . In particular one can use the body of experience gained in the study of gauge theories 
to attack the problem. With an adiabatic approximation for f(x,k), has a series expansion in powers of x. Thus 
for quasi static systems 

a M = (co + cxx + . . (159) 

The effective field strength / M „ need not always be zero. Two situations might arise: (i) the Wigner function might 
contain a logarithmic singularity, as in the case where vortices are present, and (ii) the source 7^ could contain 
components which specifically drive the macroscopic field in a given way. A simple example of the latter is the 
analogue of Rabi oscillations in the laser, in which the field oscillates between two states in a regular way. Here, this 
oscillation is driven by the source 7^ or perhaps by a pulsation of the inhomogeneity scale, and occurs from the linear 



♦7 



terms in equation (159). The current J = 2 D 4>i behaves like a magnetic influence on the system (doing no net 



work). Simplifying to the case of a (1 + 1) dimensional system, one may write 

a M ~ % = Yf^nv'tf (160) 

for constant I7I and /x, v = 0, 1. This corresponds to a harmonic 'flip-flop' motion between field and source. It is also 
directly analogous to the well known problem of Landau levels in an effective magnetic field 7|. 

The localization in spacetime resulting from the inhomogeneity scale suggests that such oscillations may take place 
locally in cellular regions. A simplified model for this is to impose periodic boundary conditions on the cells, generating 
a kind of global field coherence (this is admittedly motivated by technical simplicity rather than physical reasoning). 
One is therefore led to the study of Landau levels on the torus — a system which has been studied at some length 
|^8|,0, and will not be re-analyzed here. 
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A significant feature of the Landau problem on the torus is that the periodicity enforces a flux quantization condition 
on the field. Here this translates into the following relation: 



d Q H{k,x)diH{k,x)\~f\ = 27m (161) 

for integer n. This relation indicates that nearest neighbour cells might engage in cooperative oscillations (i.e. the 
size of cells is quantized in units of the local inhomogeneity scale). This is clearly a far less stringent condition here 
than in the case of a true periodic torus, since the inhomogeneity scale varies in space and time thus the meaning of 
strict quantization is lost. However, it indicates that one can expect a tiling of spacetime by oscillation cells. Since the 
size of the cells might be highly irregular, the tiling behaviour is most likely to be chaotic unless special geometrical 
boundary conditions can enforce a re gular ity on the field. This is an alternative expression of the behaviour deduced 
from the Green function in equation (116). 



VIII. CONCLUSION 



Schwinger's closed time path action principle has been applied to the neutral scalar meson, off-equilibrium, in the 
presence of long-range, inhomogeneous sources. The method of dispersion relations is used to find formal expressions 
for Green functions which reflect the absorbative and amplifying processes in the normal modes of the field. In the 
case of self-interacting theories, the sources can be thought of as representing <fi 2n interactions to one-loop order, 
effecting a rcsummation of the theory. The effect of rapid transport (large F M ) is to induce a change in the sign of 
the mass squared, indicating a second order phase transition and anomalous dispersion. 

If significant inhomogeneities or long range interactions exist, the field naturally forms localized cells with (to lowest 
order) a periodic relationship to the natural inhomogeneity length/time scale. This is shown from the viewpoint of the 
Green functions and by recourse to an analogy with Landau levels on the torus. Since the length scale is determined 
by non-linear considerations one can expect chaotic behaviour with islands of order (stable fixed points) along the 
approach to equilibrium. A simple analogue of Rabi oscillations in the laser is shown to arise as a leading order 
behaviour in x. 

The method used in the this work has the advantage of combining the fundamental aspects of an operator field 
theory with the usefulness of the action principle. The use of generating functional ultimately leads to functional 
integral forms, as used almost exclusively in the literature. However, the introduction of the functional integral is 
scarcely necessary using the present method and often has the undesirable effect of turning the discussion of causality 
into one of complicated paths of integration in the complex plane. 

Comparing to other works reveals both differences and similarities. Lawrie H, for example, treats the quantity 7^ 
as an explicitly written imaginary part of the spectrum of excitations. He ignores F lll but does not ignore f^. This is 
a somewhat different approximation which has a more distant relationship to classical transport theory. In fact, since 
the appearance of and in a M is identical, up to a sign, the form of dynamics might well be independent of the 
approximation used in this work — understandable as a reparameterization of an equivalent problem. Lawrie further 
considers (f> A theory and uses a renormalizatio n-lik e philosophy to determine the sources self-consistently thereby 
effecting a resummation as noted in equation (141). Calzetta and Hu Q use a variational principle to determine 
the effective action for a self- interacting boson theory. This is the same idea as in ref. ||, expressed in extremely 
aesthetic formalism and containing important insights into the subject; the solution to their method is, in practice, 
more difficult to attain however and thus results are mainly formal. Neither of these works consider the implications 
of non-local effects. Another interesting approach is the Schrodinger quantization approach in ref. pi. This makes a 
contact with the Schwinger action principle at a more subtle level and, focusing on somewhat different issues, uncovers 
features absent in other formulations of non-equilibrium physics. 

It is important to extend the present analysis to include both fermions and spin-1 bosons (true gauge fields). The 
latter is probably a difficult task in view of the problems which can arise in gauge fixing. Again, the action principle 
approach, starting from the operator field theory is likely to be the most informative approach. The appearance 
of discrete spectra and magnetic like effects makes the present work very interesting to the study of the fractional 
quantum Hall effect. In particular, the pseudo-gauge field formulation might have interesting connections with the 
statistical gauge field employed in the Chern-Simons gauge theory picture. These and other outstanding issues will 
be discussed in future work. 

I am grateful to I.D. Lawrie for helpful discussions. 
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APPENDIX A: RETARDED TV-POINT FUNCTIONS 



The retarded n-point functions are defined by 



(n = 0) R(x) = = (j>(x) 
(n = 1) R(x, xi) = = —i9(x — x\)[(j}{x), (f>(xi)] 

(n > 1) R(x, x\... x n ) = {-i) n ^ 6{x - xi)6{xi - x 2 ) . . . 0(x n -\ - x n ) x 



Pi 



[[[(j>(x), 0(xi)], <j>(x 2 )] <t>{x n )\ 



(Af) 



where Pi signifies all the permutations of the indices on X{. This is strictly only defined when all the Xi are different. 
The coincidence limit is often defined by recourse to the momentum representation. 

The retarded functions have the following properties: (i) R(x, x\ . . . x n ) vanishes if any Xi > x (with respect to the 
time), (ii) R(x, X\. . .x n ) is a symmetric function of X\ . . . x n and (hi) the retarded Green functions are always defined 
with respect to a special point x which is later than all other points. 

It may be verified explicitly that the Schwinger-Symanzik generating functional 



Z ret (x) = (TV"*) 



generates the n-point functions according to the rule 



SJ(x) 



R(x,xi ...x n ) = {-i) n — Z ret (x) 



(A2) 



(A3) 



.7=0 



The step functions are enforced by explicit cancellation of field operators for times outside the bounds of the con- 
straints. The above generating functional is clearly related to the closed time path generator, and it is easy to see 
that one may also write 



R(x,xi ...x n ) 



5 



5 



H)"^<o|o>± 



(A4) 



./±=o 



Finally, it can be observed that that the Hermite polynomials are generated by the generating functional 



H n {z) = (-l)»e- — e" 



(A5) 



and can therefore be expected to play an important role in the computation of the transformation function for a 
quadratic theory. 
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